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ABSTRACT 

The behaviour of the afterglow (AG) of gamma-ray bursts (GRBs) directly 
provides, in the cannonball (CB) model, information about the environment of 
their progenitor stars. The well observed early temporal decline of the AG of GRB 
021211 is precisely the one predicted in the presence of a progenitor's "wind" 
which resulted in a density profile oc 1/r^ around the star. The subsequent fast 
fading — which makes this GRB "quasi-dark" — is the one anticipated if, further 
away, the interstellar density is roughly constant and relatively high. The CB- 
model fit to the AG clearly shows the presence of an associated supernova akin to 
SN1998bw, and allows even for the determination of the broad-band spectrum of 
the host galaxy. GRB 990123 and GRB 021004, whose AGs were also measured 
very early, are also discussed. 

Subject headings: gamma rays: bursts 



1. Introduction 

Massive stars are observed to emit a dense stellar wind, continuously or in a sequence 
of eruptions, which presumably intensifies before their death as core-collapse supernovae 
(SNe). Indirect evidence for a large mass loss from the progenitor stars during their late- 
time evolution is provided by the absence of H lines (and He lines) in the spectra of SNe of 
Type lb (Ic). In Type II SNe, direct evidence for a wind emitted shortly before the explosion 
is provided by the subsequent observations of emission by the circumburst material of narrow 
optical lines (e.g. Salamanca et al. 1998; 2002; Fassia et al. 2001), radio waves and X rays (for 
a recent review, see Chevalier 2003). All these observations indicate that the circumstellar 
density profile is p^j oc r~^, as expected from a quasi-steady stellar wind. 

Gamma ray bursts (GRBs), in the cannonball (CB) model, are associated with the death 
of massive stars in core collapse SN explosions (Dar & De Rujula 2000a,b; Dado, Dar & De 
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Riijula 2002a,b). Clear indications for this association are the effect of the progenitor's wind 
on their early-time afterglow (AG) and the presence of a supernova in their late-time AG. In 
a CB-model analysis, there is very convincing evidence of a GRB-SN association (Dado et 
al. 20002a-e; 2003). The CB-modcl fingerprints of the massive wind of the progenitor star 
have also been spotted in the AGs, but never before with unquestionable certainty. 

In this paper we fit the broad-band optical afterglow of GRB 021211 in the CB model, 
and extract information about the circumburst density profile, the associated SN and the 
host galaxy. We also rediscuss GRB 990123 and GRB 021004. In these three cases the AG 
was observed in the optical band early enough after the onset of the burst for the CBs to be 
still moving through the progenitor's wind. The observed early AG decline, roughly t~^'^, 
is the CB-model expectation^ for a wind density profile oc r~^. In the well measured early 
AG of GRB 021211, the agreement between the observations and the CB-model prediction 
is extremely satisfactory. 

In the CB model we assumed GRBs to be associated with SNe similar to SN1998bw 
(Dar and De Riijula, 2000a). In subsequent work (Dado et al. 2002a-e; 2003) we found 
that, indeed, in all cases of GRBs of known z in which such a SN was in practice visible (all 
cases with z < 1.12), the CB-model fit to the corresponding AG disclosed its presence. GRB 

021211, at ^ = 1.006, is one more example: the fit clearly requires such a SN contribution, 
tentatively reported by Fruchtcr et al. (2002) and Testa et al. (2003). The CB model 
has so few parameters that we could determine via our fits the unknown contribution of 
the host galaxy of GRB 021211 to the different optical bands. The resulting broad-band 
spectrum snugly resembles that of other GRB host galaxies (e.g. Gorosabel et al. 2003) and 
star- forming galaxies at a similar z (Pruchter et al. 2002). 

In the CB model "dark" GRBs without AGs are not expected. But the AGs' decline is 
expected to be often fairly fast, in particular if the interstellar medium (ISM) density is rela- 
tively high, making their detection difficult or, in the bygone days of slow GRB localization, 
nearly impossible; see Fig. 6 of Dado et al. (2002a) and its discussion. With faster detectors, 
more of these faint, fast-declining "quasi-dark" GRBs ought to be found, and GRB 021211 
is an example. The detailed shape of the fast AG dechne observed from 20 to 150 minutes, 
roughly as 1/t, is well fit with an ISM whose density, beyond the reaches of the wind, is 
relatively high and approximately constant. 



■^We have claimed before that this behaviour ought to be ^ t~^. In the case of GRB 990123, discussed 
in Dado et aL 2002a, that was due to our then-incomplete understanding of the AGs' broad-band spectra, 
which we improved in Dado et al. 2002b. For GRB 021004 it was an error, which we correct here, although 
it is not relevant in practice (the early data on this AG are not sufficient to distinguish the two declines). 
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2. GRB 021211 

The gamma ray burst 021211 was detected on December 11.471227 UT, 2002, with the 
HETE FREGATE, WXM, and SXC instruments. It was a relatively long-duration (> 5.7 s 
in the 8-40 keV band), single-pulse GRB with a peak flux of 8 x 10^^ erg cm^^ s^^ and a 
fluence of ~ 10~^ erg cm~^ (Crew et al. 2002). Its early-time optical AG was detected by 
the robotic telescopes RAPTOR (Wozniak et al. 2002), KAIT (Li et al. 2003) and S-LOTIS 
(Park et al. 2002), 90, 108 and 143 s after the burst, respectively. A fit to the KAIT first 9 
data points, in the interval of 2.2-6.5 minutes after burst, showed a temporal decay with a 
power-law index of —1.60 ± 0.02 and a fit to later data, taken 20 to 150 minutes after the 
GRB, gave a slope of —0.96 ±0.04 (Chornock et al. 2002), consistent with the decay seen by 
Price and Fox (2002a,b) who were first to report the detection and the precise localization 
of the optical AG. Follow-up observations in the optical, NIR and radio bands were reported 
by Fox et al. (2003) and by other groups in the GCN. The redshift of the host Galaxy, 
z — 1.006 ± 0.002, was measured with VLT by Vreeswijk et al. (2002) and confirmed by 
Delia Valle et al. (2003). A fiattening of the fight curve due to the host galaxy and a possible 
contribution from an associated SN were observed with HST 7 and 14 days after burst in 
the BVIH bands (Fruchter et al. 2002) and with VLT-UT4 between days 30 and 35 in the R 
band (Testa et al. 2003). Well-sampled optical AG observations have been recently reported 
by Li et al. (2003). The published observational data on the broad-band optical AG of GRB 
021211 are shown in Figs. 1 and 2. 



3. The CB model 

In the CB model bipolar jets of CBs are launched axially in core-collapse SNe, with 
initial Lorentz factors 70 = (9(10^). The CBs are assumed to be produced in an unstable 
accretion process, as in quasars and micro-quasars and, as in SS 433, to be made of ordinary 
matter^. Crossing the SN shell and the progenitor's wind with a large 7, the front surface of 
a CB is coUisionally heated to keV temperatures. The quasi-thermal radiation it emits, when 
no longer absorbed by the intervening matter, and boosted and coUimated by its relativistic 
motion, is a single 7-ray pulse in a GRB. The cadence of pulses refiects the chaotic accretion 
and is not predictable, but the individual-pulse temporal and spectral properties are (Dar 
and De Riijula, 2000b; for recent reviews, see De Riijula 2002; Dar 2003). The ejected CBs, 
as observed in //-quasars, are assumed to contain a tangled magnetic field. As they plough 



''Balmcr H and He lines (e.g., Eikenberry, et al. 2001) and the Ka line of Fe (Migliari et al. 2002) were 
detected in the mildly relativistic CBs of this /U-quasar. 
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through matter, they gather and magnetically scatter its constituent protons. The re-emitted 
protons exert an inward pressure on the CBs, which counters their expansion and makes them 
reach an asymptotic radius i?^^, in minutes of observer's time (Dado et al. 2002a). The 
electrons swept in by the CB in its voyage through the wind and ISM are Fermi-accelerated 
in the enclosed magnetic maze and cooled by synchrotron radiation. Shortly after the GRB, 
the optical radiation from the CB is dominated by synchrotron emission from these electrons. 
So far, the CB model was very successful in fitting the observed broad-band AGs of all GRBs 
of known redshift (Dado et al. 2002a-e, 2003). 



4. The fate of a cannonball 

Let Up be the baryonic number density of the circumburst material or the ISM, both 
dominated by protons. A spherical CB of radius R^^ flying through this material sees, in 
its rest system, an incoming fiux of protons entering it with a total momentum per unit 
time n ~ n^mp 7^ Trit!^^, with 7 = ^{t) the diminishing CB's Lorentz factor. If the CB's 
enclosed magnetic field randomizes these protons so that they are re-emitted isotropically, 
the CB's surface is subject to an inward pressure P = 11/(4:7: R"^^), which is independent 
of R^g- In Dado et al. 2002a we have assumed that this pressure is sustained by the 
equal and opposite pressure 5^/(8 tt) of the CB's enclosed magnetic field (which is thereby 
determined) resulting in an equilibrium situation for a CB of approximately constant R^^g- 
We also assumed that R^g , for a CB originally expanding at a transverse velocity Po c, 
can be estimated as the maximum expansion radius attained as the pressure P opposes the 
expansion, that is i?^^ ~ [3 Ai'^g/5Q/(2 TrUp 7q)]^/^, with N^^ the CB's baryon number. 

To test the above bold assumptions, we have fitted the model to all broad-band AGs of 
GRBs of known z, with satisfactory results. In Dado et al. (2002a) we have also tried a model 
with continuously expanding CBs, and found it to be completely inadequate, supporting the 
ansatz of CBs with constant radius. The estimated value of S ~ 10 Gauss (for Hp — 10~^ 
cm~^ and 70 = 10^) is adequate for the description of the data, but our original guess 
R^g — 0(10^^) cm turned out to be an overestimate by at least one order of magnitude 
(Dado et al. 2002b). In this paper we subject the CB's defining property — that they reach 
an approximately constant radius — to a severe test, by studying in detail its consequences 
for very early AGs. 

Let 6 be the angle between the direction of a jet of CBs and the observer. The Doppler 
factor of the light the CB emits is well approximated by S{t) ~ 2^[t)/{l + 9'^^{ty) in the 
domain of interest for GRBs: large 7 and small 9 (typical values ensuing from our fits are 
Jo 1/9 10^). The relation between time — as measured by the observer — and distance 
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— as travelled by the CB — is: 

dx^^^cdt. (1) 
The ambient protons that a CB scatters in the interval dx slow it down by an amount^: 

i,= -^%M!d., (2) 

where N^^g is the CB's baryon number, for which our reference value is 6 x 10^°. The function 
7(t) can be explicitly found by quadrature in the two cases of interest here: a density profile 
n{x) — nyj{xw/xY (for a "windy" neighbourhood) and a constant density (an adequate 
approximation as the CBs get further away into the ISM). 

For an ISM of constant density, ^{t) depends on d, on the initial 7 = 70 as a CB exits 
the denser wind domain, and on Xoo — N^j^/{Trnp R'^^), a deceleration parameter. The value 
of j{t) is the real root of the cubic: 



T 7o 



1 1 

,7 7o, 



(3) 



(1 + z) X 



It takes a distance xi/2 — Xoo/jo for ^{t) to descend to 7o/2. Our fitted values of 2:1/2 are in 
the range 0.1 to 1 kpc, corresponding to AGs that fade in the observed characteristic times 
of order (1 + z) Xi/2/(cSo) ~ days. This range may be affected by an observational bias: 
AGs with a smaller Xoo (e.g. those with CBs moving in a relatively high density ISM) decay 
faster in time and are "quasi-dark": harder to detect. GRB 021211 is one such case. 

We do not report the function ^{t) for the small values of t corresponding to the CBs 

crossing the parent stellar wind, since, for typical parameters, the fractional energy loss 
A7/7 is negligible in that interval. The "canonical" stellar wind of a very massive star has 
a rate M^, ~ 1O~^M0 y~^ and a velocity f„, ~ 100 km/s, so that p{x) ^ {x^u)/xY with 
Pwx1j^M,uj/ {4:11 v.^^) ^'0.1 X 10^'^ g cm^^, or n„^^,T,^~3 x 10^'' cm~^. The actual observations 
of Myj and Vy, span an order of magnitude around these canonical values. It follows from 
Eq. (2) that the condition A7/7 ^ 1, for a wind profile extending from to a;/ 3> Xj, is 
-^CB ^-^cB^i/(^"'w^^7)- ^'^^ observations starting at = 100 s after burst, z — 1 and 
typical ^ — 5 — 10^, Xi — C7(5 Ai/(1 -\- z) — 0.5 pc, and the constraint is R^g < IQi^^ cm, 
which we have found to be amply satisfied (Dado et al. 2002b). 

The fact that A7/7 is typically small as the wind is crossed has, as we shall see anon, 
an important consequence: the shape of the early optical AG locally and directly reflects the 
shape of the circumburst density profile. 



^We call X and not r the distance from the progenitor, not to insinuate a hypothesis of spherical symmetry. 
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5. The GRB afterglow in the CB model 

The AG — the persistent radiation in the direction of an observed GRB — has three 
origins: the ejected CBs, the concomitant SN explosion, and the host galaxy (HG). These 
components are usually unresolved in the observed "GRB afterglows", so that the light 
curves and spectra are measures of the corresponding cumulative energy flux density: F^^ = 
^cBs + ^SN + Fhg- observed cases, but one (GRB 021004, discussed in Dado et al. 

2003), it is sufficient to approximate the ensemble of CBs in the AG phase as a single (or 
dominant) CB. 

The contribution Fsn of the SNe is approximated by that of SN1998bw (Galama et al. 
1998), displaced to the GRB's redshift^ In the CB model the pair GRB 980425 & SN1998bw 
is in no way exceptional (but for the accidentally small z and large 9 of the GRB). Thus, 
the use of SN1998bw as a candidate GRB-associated standard candle makes sense. 

The optical AG of a CB is dominated by synchrotron radiation from the ISM electrons 
that penetrate in it. We argued in Dado et al. (2002b) that these electrons are Fermi- 
accelerated in the CB-enclosed magnetic maze and cooled by synchrotron radiation to a 
distribution (in the CB's rest system) due/dE oc below the injection bend energy at 
which they enter the CB (Eh = ^(t) rrieC^); and to a distribution dn(,/dE oc above 
-Ef,, with p = 2 in analytical approximations and p = 2.2 in numerical simulations. The 
emitted synchrotron radiation has a corresponding double power-law form vdn^/dv oc 
with a power-law index ai ~ 0.5 well below a bend frequency Ub, and ah = p/2 ^ 1.1 
well above it. Broad-band fits to the data on the AGs of all GRBs of known z result in 
ai 0.6 ± 0.1, afi ~ 1.1 ± 0.1. In practice a; is less well determined than ah, the CB-model 
fits being good even for a fixed ah — 1-1 ■ We attribute this to the fact that the extracted ah 
is sensitive to the AG light-curve at relatively late times, when CBs are typically hundreds of 
parsecs away from the progenitor and the absorption in the host, which is hard to ascertain, 
may be minimal. 

The value of Ej, and the estimate of the CB's magnetic field described in the previous 
section allow us to estimate the bend frequency of the synchrotron-radiation spectrum: 



1.87 X 10^j{t)f6{t) 

iT~z 



np{x) 



1/2 



Hz. 



(4) 



In the observer frame, this radiation is Doppler-boosted and coUimated by the relativistic 
motion of the CB, and redshifted by the cosmic expansion. Its explicit form, extending 



^The cosmological parameters we use are: Hq = 65 km/(s Mpc), f^M = 0.3 and = 0.7. 
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from radio to X-ray frequencies, is given in Eqs. (4,6,8) of Dado et al. 2002b; here we only 
report its behaviour near the optical domain, which is simple above and below the injection 
bend. Let a = («;, ah) = (0.5,1.1) be the predicted power indices at z/ <C z^fe and u^Uh, 
respectively. The time and frequency dependence of the energy fluence in these limits is: 

oc A'(u, t) [7(i)]3"-i [S{t)]^+'' i/-" , (5) 

where rie = ng{x) np{x) is the electron density along the CBs' trajectory and A'(z/, t) 
corrects for absorption in the host galaxy and in ours, its possible time dependence originating 
in the kiloparsec length of the CBs' trajectory in the host galaxy. Self- absorption in the CB 
is irrelevant at the optical and NIR wavelengths relevant here. In the CB model the early 
AG is dominated by the last significant pulse (or CB) of the GRB, so that t in Eq. (5) is the 
time after that pulse. The nuance may be significant at very small t. 

In deriving Eq. (5) we have assumed that a fixed fraction of the energy-deposition rate 
by ISM electrons in a CB (vri?^^ nemgC^T^, in the CB's system) is re-emitted as the AG. 
The AG spectrum is a function of u and Ub, so that at fixed u it depends on rig = rip via u^, 
as in Eq. (4). Thus a result which is not linear in rie, and peculiar at first sight. 

In our fits we assume as the circumburst density profile a constant plus a "windy" term: 
n{x) = no + nw{xw/xY = no [1 + [x/xf]. (6) 



6. The very early optical AG 

For typical parameters, the first few parsecs of a CB's voyage — as it crosses the "wind" 
material emitted by the progenitor star — are seen by an observer in the first few minutes of 
the AG. During that "early" time, the ambient density is high enough for z/6, as in Eq. (4), 
to be comfortably above the optical frequencies, so that a = a; = 0.5 in Eq. (5). We have 
seen that, at "early" times, a CB has its Lorentz factor insignificantly changed by collisions 
with the wind material, so that the observer's time and the CB's travelled distance, as in 
Eq. (1), are strictly proportional. Thus Eqs. (5,6) collapse to: 

F,^n^^+-y^o,[l + {t/tf]"\ (7) 

Had we used the observed a; = 0.6 ±0.1, as opposed to the naive theoretical a; = 0.5, 
we would have obtained F^, ~ t~^-^ at veary early times, in close agreement with the fit by 
Chornock et al. (2002) for GRB 021211, and the observations of Akerlof et al. (1999) for 
GRB 990123 and Fox et al. (2002) for GRB 021004. But we shall see that Eq. (7) provides 
an excellent fit to the data. The observed deviations relative to the smooth curve of Eq. (7) 
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reflect the variations of the wind that generated the profile. In the case of GRB 021211, 
these variations are at the ~ 10% level (Li et al. 2003). 

7. The overall AG of GRB 021211 in the CB model 

The copious early data on GRB 021211 obtained by Li et al. (2003), along with all 
other data communicated in the GCN notices quoted in section 2 and in Fox et al. 2003 are 
shown in Fig. 1 for the R-band, and in Fig. 2 for the broad-band optical (IRVB) and NIR 

(HKJ) passbands. In fitting these results we neglect the unknown extinction in the host, due 
to the lack of relevant spectral information^. For this GRB the R-band data are relatively 
so copious that it is useful to perform separate R-band and broad-band fits. 

In our fit, the density is that of Eq. (6) and the limiting exponents in the synchrotron 
spectrum oc are set to their theoretical values {ai,ah) = (0.5,1.1). The fit is made 
with the complete analytical broad-band predictions Eqs. (4,6,8) of Dado et al. 2002b, of 
which Eq. (5) are the limiting behaviours and Eq. (3) is the deceleration law. The best fitted 
parameters for the R-band fit shown in Fig. 1 arc 70 = 262, 9 = 1.76 mrad (i.e. 6q = 431) 
and Xoo = 4.7 kpc, with 51 d.o.f. and x^/(d.o.f.) = 0.97 if the out-lying point at t~2 days 
is eliminated. The constant contribution of the host galaxy was left as a free parameter, for 
which the fit returned F^^ — 0.34 //Jy. The contribution of a standard-candle SN1998bw at 
the GRB position is discernible in Fig. 1, in which F^^ has been subtracted. 

We have also tried a fit with an arbitrary power a in the windy term: n^^x^/ xY , 
resulting in a best fitted a — 1.92 ± 0.03, close to the "canonical" expectation a — 2. Thus, 
the agreement of the fits with the canonical "windy" power is not a consequence of the 
scarcity of data. Also, at early times, the fit is totally insensitive to the other parameters 
(70, 6 and Xoo)- Thus, the result for a is robust: the signature of a wind is there. 

The overall normalization of the AG is also a free parameter in our fits; it is proportional 
to the ambient density n{x), to the CBs' cross sections, to their number, and to the unknown 
fraction of the energy of the gathered ISM electrons that is re-emitted as the AG. Thus, 
extracting a value of n{x) from the fit normalization requires extra assumptions. Yet, the 
value of n{x) at some reference distance (e.g. tiq, its limiting value at large distance) can be 
extracted from the fact that we use it as a free parameter in the expression of Eq. (4) for 
V},. The fitted result is Uq = (2.97 ± 0.22) cm~^. Combined with the fitted value x = 1.2 



^An unabsorbed SN1998bw-like contribution fits the late-time broad-band data, indicating that host 
extinction is not large. 
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pc in Eq.(6), this yields py, = mpUy, x"^ — (6.8 ± 0.5) x 10^^ g cm ^, compatible with the 
expectation, quoted in section 4, for a "canonical wind": 5.1 x 10^^ g cm~^. 

The wide-band fit is shown in Fig. 2; it has essentially the same fit parameters as the R- 
band fit. The contributions F^^ of the host galaxy in the different bands, resulting from the 
fit, are 0.50, 0.34, 0.23 and 0.23 //Jy for the IRV and B bands, respectively. This extracted 
spectrum, and in particular its flattening above the V frequencies, is compatible with that 
of other GRB hosts (e.g. Gorosabel et al. 2003) and star-forming galaxies at z ~ 1 (e.g. 
Fruchter et al. 2002). 

The comparison we have described between the CB-model expectations and the data for 
GRB 021211 is simply spectacular. The fits are good, particularly that to the shape of the 
R-band light curves. The agreement between the expected and extracted ambient density 
profiles is eery. A word of caution, however, is necessary. The number of parameters in 
our fits is high: the "ambient" parameters describing the host galaxy (F^^) and the density 
profile (no and x), the observer's viewing angle 9, and the 3 quantities that are specific to 
the CBs: the overall AG normalization, 70 and Xoq. Even if the fitted ambient parameters 
turn out to have the expected values, a model with this many parameters has to be very 
wrong not to describe a very simple surface: the fiuence as a function of frequency and time^. 
Moreover, the data are not cross-calibrated, their systematic errors are unknown, the model 
is no doubt a simplification of a very complicated phenomenon, and the parameters are not 
uncorrelated (except at the early times of particular interest here). As a consequence of all 
this, the "formal errors" of the parameters (as given by the standard program MINUIT) are 
no doubt underestimated. Whence our habit of not always making the errors explicit. 



8. The radio AG of GRB 021211 

In the CB model the radio AG at early time is suppressed because of limb darkening 
and the finite time it takes the electrons gathered by a CB to cool radiatively and reach 
an energy at which they emit synchrotron radiation of a given typical frequency (Dado et 
al. 2002b). At later times the radio emission, like the optical and X-ray emission, becomes 
proportional to a high power of the Lorentz factor. In the case of the radio AG of GRB 
021211, only upper limits (or perhaps marginal detections) have been reported, by Hoge et 
al. (2002) and Fox et al. (2002). The later authors, who analyze the AG in detail, consider 
the absence of a stronger radio signal to be a problem. In the CB model, it is not. 



It is possible to construct models that disagree with observation or are unphysical, but have even more 
parameters. For the case of GRB 991208, see section 2 of Dado et al. 2002e. 
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The wide-band spectrum of AGs is, in the CB model, very simple. Its only feature 
— besides the bend frequency that we have discussed — is due to self-absorption in the CB, 
and it is described by a single parameter z/^ (Dado et al. 2002b). In the case of GRB 021211 
we cannot fix this parameter, having no single secure radio signal. The best we can do is to 
show, by way of example, what the radio signal would be if the value of Ua was similar to 
the average for other radio AGs - 1 GHz) or to the value for GRB 021004 (0.98 GHz), 
which is otherwise akin to GRB 021211, and is very well described by the CB model (Dado 
et al. 2003). The result of this exercise is shown in Fig. 4, and it is compatible with the 
observational limits or marginal detections: the largest "signals" reported in Fox et al. 2002 
are 45 ± 23 and 60 ± 38 /xJy, at i ~5, 10 days, both at 8.46 GHz, while Hoge et al. (2002) 
obtain a 3a upper limit of 7.5 mJy at 347 GHz, on day t~l. 



9. The GRB proper 

In the CB model, the fluence F from a GRB viewed at a small 9, is amphfied by a huge 
factor 5q, due to Doppler boosting and relativistic coUimation (Dar & De Riijula 2000): 

_ (1 + z) 51 

where is the total energy in photons emitted by CBs in their rest system. The total 
"equivalent spherical", or would-be isotropic energy, E'"", inferred from the observed fiucnce, 
is a factor (5o)'^ larger than E^. In Dado ct al. 2002a we deduced that the E^ values of the 
GRBs of known z span the surprisingly narrow^ interval IQ^^^''-^ erg, the spread in E'"'" being 
mainly due to the spread in their values of 6o (deduced from the fits to their AGs). For GRB 
021211, the CB-model expectation is ^ S^E^ 0.8 x IQ^'^^'^-^ erg, in agreement with 
the observed E^^° f« 1 x 10^^ erg, deduced from its measured redshift {z — 1.006; Vreeswijk 
et al. 2002) and fiuence in the 8-40 keV band (~ 10"^ erg cm-^; Crew et al. 2002). 



10. GRBs 990123 and 021004 

The AG of these two GRBs was also observed particularly early, by Akerlof et al. (1999) 
in the case of GRB 990123. We have recalibrated their data assuming V — R — 0.28 mag 
(our predicted ai — 0.5), rather than a constant colour. We have also recalibrated the later- 
time data compiled in Castro-Tirado et al. 1999; Kulkarni et al. 1999; Galama et al. 1999 



^GRBs in the CB model are much better standard candles than in the standard model (Frail et al. 2001). 
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and Pruchter et al. 1999, assuming R — r — 0.015 mag {ai — 0.5). The R-band results of 
our broad-band fit are shown in Fig. 4, in which t = corresponds to the start of the last 
(second) prominent pulse in the GRB, 37s after trigger (Briggs et al. 1999). There is an 
excellent agreement between the early data tracing the circumburst windy density and the 
prediction of Eq. (6). A fit with an arbitrary power a in the windy term, n„,(x„,/x)", results 
in a best-fitted a = 1.98 ±0.06, again in agreement with the "canonical" expectation: a = 2. 

Though the data have been recalibrated, the parameters we fit to the AG of GRB 990123 
are close to those in Dado et al. 2002a,b. In particular 70 = 1204 and 60 = 1630 — determined 
from the shape of the AG — are the highest for any GRB of known z, implying that the 
magnitude of the GRB and AG fiuences should also be record-breaking, as they are (even 
for this GRB, the value of Ej in Eq. (8) is not exceptional). Since 7 and S are so large, the 
bend frequency uj, of Eq. (4) "crosses" the optical band very late in time, where the data are 
quite imprecise. As a consequence, the R-band and wide-band fits are very insensitive to uj,, 
and we cannot reliably use its value to determine the absolute magnitude of the ambient and 
wind densities. If we assume that the density profile prevails until Uw declines to 

the typical superbubble or halo value, Up ~ 10~^, we find that pw x"^ = mpfipxl, ~ 4.2 x 10^^ 
g cm~^, again in agreement with the expectation, quoted in section 4, for a "canonical wind" : 
5.1 X 10^^ g cm~^. The windy profile would in this case extend all the way to x^, ~ 50 pc, 
yielding a total mass loss of ~ 40 Mq during the fife of the very massive progenitor star. 
Such a large total mass loss is not unusual for the progenitors of SNe of Type lb and Ic. 

The CB-modcl fit to the R-band AG of GRB 021004 is shown in Fig. 5. The data arc 
those measured and compiled by Fox et al. 2002; Holland et al. 2002; Pandey et al. 2002 
and Bersier et al. 2003. Once more, the agreement between the prediction of Eq. (6) and 
the early data is good, though these data are neither as abundant nor as close to the GRB 
trigger as for the other two cases we discussed, implying that we cannot claim to have clearly 
spotted, for this GRB, a typical windy profile. The fitted result for the asymptotic density 
is no = (4.08 ± 0.17) x 10~^ cm^'^. Combined with the fitted value a; ~ 5.1 pc in Eq. (6), 
this yields Pwxij = mpnwxlj ~ 1.7 x 10^^ g cm~^, a bit over an order of magnitude below 
the central "canonical" expectation. 



11. Discussion and conclusions 

It is, as usual, interesting to compare the CB model with the standard fireball or blast- 
wave models of GRBs. For GRB 021211 a detailed SM analysis is given in Fox et al. 2003. 
The relevant issues at hand concern the signatures for "windy" neighbourhoods, and the 
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predictive power for early AGs^°. 

In the standard model (SM) the scarcity of AGs indicating a windy circumburst density 
distribution is a problem, admitted even by its staunchest defenders (Piran 2001; Price et 
al. 2002). In the CB model all the observed AGs that should show a windy signature do. 
Indeed, the relation Eq. (1) between observer's time and CBs travelled distance imphes that, 
for the typical 7o ~ (^o ~ 10^, the AG of a GRB at a cosmological distance must be "caught" 
within the first few minutes for the CBs to be still travelling close enough to the parent star, 
so that the contribution of its wind to the ambient density is still significant and observable. 
Only the three GRBs that we have discussed pertain to this category. 

In the SM the fast-decreasing "early" optical AG, as first observed in GRB 990123, is 
generally attributed to a reverse shock (Meszaros and Rees 1999). For GRB 021211, this 
is the interpretation spoused by the observers (Fox et al. 2003; Li et al. 2003). In this 
view the index of the approximate power law of the decay of the early AG is not fixed, 
and its normahzation is unrelated to that of the late AG. In the GB model, contrariwise, 
the early optical AG does not have its own separate origin: it is made by the very same 
mechanism as the late AG. The temporal shape of the AG is not a succession of power laws 
with breaks between different mechanisms or regimes. This temporal shape is not arbitrary: 
at early times it is directly related to the shape of the circumburst density profile, as in 
Eq. (7), and this relation is fully vindicated by the data for the expected windy profile. 
The AGs shown in our figures depend, in the "early" domain, on only one parameter: the 
absolute normalization. The ambient density that defines the early temporal shape can also 
be extracted from the CB-model fit, via the bend frequency of Eq. (4), and its magnitude 
turns out to coincide with the expected density for the canonical wind of a heavy star. 

Fast-declining optical light curves — masquerading as "dark" AGs — are expected in the 
CB model (Dado et al. 2002a), and GRB 021211 is an example. The CB-model fit to its 
AG has allowed us to infer the presence of an associated SN akin to SN1998bw, as for all 
other GRBs where the SN could in practice be detected. The fit also yields the "colour" of 
the host galaxy, which is the expected one. 

We have seen how, concerning the early optical AGs, the CB model is very predictive 
and successful. In previous works we have shown this to be the case also for all properties of 
the AGs of all GRBs of known redshift: they can all be described in a simple, analytical and 
parameter-thrifty manner, which does not involve multiple choices, multiple mechanisms and 
multiple exceptions. To put it in a way with which even the most allegiant apologists of the 



-^'^For Fox ct al. (2003) the non-observation of a radio AG is also significant, suggesting that "the burst 
may have suffered substantial radiative corrections" . No such specific cure is necessary in the CB model. 
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SM might agree: the contrast between the CB model and the standard model is striking. 
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Fig. 1. — The optical observations of the R-band AG of GRB 021211, and their CB-model 
fit. The ISM density is a constant plus a "wind" contribution decreasing as the inverse 
square of the distance. The two contributions are equal at x ^ 1.2 pc, a distance reached 
by the CBs in an observer's time t ~ 0.025 days after burst. The data are those reported to 
date, in the GCN notices quoted in section 2, in Fox et al. 2003 and in Li et al. 2003. The 
contribution of a SN1998bw-like SN at the GRB position is discernible at late times. The 
host galaxy's contribution, which was fitted, is subtracted in this plot. 
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Fig. 2. — The wide-band CB- model fit to the optical and NIR observations for GRB 021211. 
For clarity, only the better-sampled I, R, V and B bands are shown, scaled by 10, 1, 1/10 
and 1/100, respectively. The contribution of a SN1998bw-like SN at the GRB position is 
discernible at late times. The host galaxy's contributions, which were fitted, are subtracted 
in this plot. 
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Fig. 3. — The optical R-band AG of GRB 021211, superimposed on the predictions for the 
radio AG at the various labelled frequencies, for an assumed t'a = 1 GHz. 
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Fig. 4. — The optical observations of the R-band AG of GRB 990123, and their CB-model 
fit. The ISM density is a constant plus a "wind" contribution decreasing as the inverse 
square of the distance. The contribution of the host galaxy has been subtracted. The small 
extinction in the Galaxy (Schlegel et al. 1998) was neglected. 



- 20 - 



10 



10 



10 



10 



_ . • 1 1 1 1 1 M 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 M 1 1 


1 1 1 1 1 M 1 1 1 1 1 1 1 M L 


TOTAL 




'■■■ii. 



\ 


R-BAND 






- CB2 \_ 

- 


- 

\ 


CB1 


\ 

. \ - 


GRB 021004 


i. ^ 
••. i '-- 

t 

' I 


1 1 1 


■ t 

■ 1 

1 1 1 1 1 'i 1 1 1 



10 



10 



10 1 

DAYS AFTER BURST 



10 



Fig. 5. — The optical observations of the R-band AG of GRB 021004, and their CB-model 
fit with two CBs, whose individual contributions are depicted along with the total (Dado 
et al. 2002e). The ISM density is a constant plus a "wind" contribution decreasing as the 
inverse square of the distance. These two contributions are equal at x ^ 3 pc, a distance 
reached by the CBs in an observer's time t ~ 0.01 day after burst. The data are the same as 
the ones quoted in Dado et al. 2002e. The contribution of a SN1998bw-like SN at the GRB 
position and the host galaxy's contribution were subtracted in this plot. 



